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Abstract. We derive integral formulas, involving the Airy function, 
for moments of the time a two-sided Brownian motion with parabolic 
drift attains its maximum. 



1. Introduction 



Let W(t) be a two-sided Brownian motion with W(0) = 0; i.e., (W(t)) t >o 
and (W(— t))t>0 are two independent standard Brownian motions. Fix 7 > 
0, and consider the Brownian motion with parabolic drift 



We are interested in the maximum 

M 7 := max WJt) (1.2) 

' -oo<i<oo ' 

of W 7 (which is a.s. finite), and, in particular, the location of the maximum, 
which we denote by 



attained at a unique point, so is well-defined a.s.) 

The parameter 7 is just a scale parameter, see Section [21 so it can be fixed 
arbitrarily without loss of generality. 

The distribution of was called Chernoff 's distribution by Groeneboom 
and Wellner (l3 | since it apparently first appeared in Chernoff 0. It has 
been studied by several authors; in particular, Groeneboom [id. 1 1 11 ] gave a 
description of the distribution and Groeneboom and Wellner [141 ] give more 
explicit analytical and numerical formulas; see also Daniels and Skyrme 0]. 
It has many applications in statistics, see for example Groeneboom and 
Wellner [l4| and the references given there, or, for a more recent example, 
Anevski and Soulier 0|. 

The purpose of the present paper is to give formulas for the moments 
of V-y in terms of integrals involving the Airy function Ai(x). (Recall that 
Ai(x) satisfies Ai"(x) = xAi(x) and is up to a constant factor the unique 
solution that tends to as x / +00. See further [l|, 10.4].) 

Date: 18 September 2012. 

2010 Mathematics Subject Classification. 60J65. 

Partly supported by the Knut and Alice Wallenberg Foundation. 



W 7 (t) := W{t) --ft 2 . 



(1.1) 



in other words, = t 
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All odd moments of ~V~ vanish by symmetry, and our main result is the 
following formula for the even moments, proved in Section [5j (The special 
case n = 2 is given by Groeneboom 

Theorem 1.1. For every even positive integer n, there is a polynomial p n 
of degree at most n/2 such that 

In particular, the variance of Vj is 

_ _ 2 -V3 7 -V3 H- _ 2 -^.y, „ 

6vri 7_ ioo Ai(z) 2 6vri 7^ Ai(iy) 2 V ; 

The integrals are rapidly converging and easily computed numerically by 
standard software. 

The polynomials p n (z) can be found explicitly for any given n by the 
procedure in Section [5j but we do not know any general formula. They are 
given for small n in Table HJ See further the conjectures and problems in 
Section El 

Numerical values of the first ten absolute moments are given by Groene- 
boom and Wellner [3]; the first four were computed by Groeneboom and 
Sommeijer (1984, unpublished). (Their values for the even moments agree 
with our formula.) 

The maximum M 7 also appears in many applications. Its distribution is 



given in Groeneboom lid . I 111 . Il2l j and Daniels and Skyrme |2j , see also, for 
example, Barbour 5|]7 Daniels 0,0], Janson, Louchard and Martin-L6f 



[151 ] . (Groeneboom [ll| describes even the joint distribution of the maximum 
M 7 and its location V~, see also Daniels and Skyrme [9].) Formulas for the 
mean are given by Daniels and Skyrme 0], see also Janson, Louchard and 
Martin-L6f [lj|; in particular 

2 -2/3 7 -i/3 

2vri 7_ ioo Ai(» 
Comparing (|1.6p and (jl.5p . we find the simple relation [1J 



EM 7 = 7^: / a:/.A9 dZ - ( L6 ) 



EV 2 = — EM r (1.7) 
Since M 7 = Wy^ — 7V?, this implies that, at the maximum point, 

EWy 7 = |EM 7 and E 7 F 7 2 = i EM 7 = J E Wy 7 . (1.8) 

It would be interesting to find a direct proof of these simple relations. 

Formulas for second and higher moments of M 7 are given in Janson, 
Louchard and Martin-L6f [151 ]; however, they are more complicated and do 
not correspond to the formulas for moments of \~ in the present paper. It 
would be interesting to find relations between higher moments of M 7 and 
moments of Ky. 
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Let us finally mention that the random variable is the value at a fixed 
time (0, to be precise) of the stationary stochastic process 

Yy(x) := argmaxj {W(t)- 7 (t-x) 2 ), (1.9) 

which is studied by Groeneboom [10 , [HI, [HI] . It would be interesting to find 
formulas for joint moments of V^{x), in particular for the covariances. 

2. Scaling 

For any a > 0, W(at) = a 1//2 iy(i) (as processes on (—00,00)), and thus 
Vj = a argmax(iy(ai) — 'j(at) 2 ) 

= aargmax(a 1 / 2 VF(t) - a 2 jt 2 ) = aV a s/2 r (2.1) 

The parameter 7 is thus just a scale parameter, and it suffices to consider 
a single choice of 7. Although the choices 7 = 1 and 7 = 1/2 seem most 
natural, we will use 7 = 1/^/2 which gives simpler formulas. We thus define 
V := an d have, for any 7, 

Vy 4 2 - 1 /3 7 -2/3y_ (2 2) 

Similarly, 

M 7 = 7i 1/3 7- 1/3 M 7 i ( 2 -3) 

for any 7, 71 > 0. 

3. First formulas for moments 



By Groeneboom [111 . Corollary 3.3], V has the density 

f{x) = ig(x)g(-x), (3.1) 

where g has the Fourier transform, see [ill . (3.8)] (this is where our choice 
7 = 2 -1 / 2 is convenient), 

/oo 2^/ 2 
e ltx g{x) dx = ■— y -00 < t < 00. (3.2) 

Note that |Ai(it)| -)• 00 rapidly as t -)• ±00, see 0, 10.4.59] or [H, (A. 3)], 
so g(t) is rapidly decreasing. In fact, it follows from the precise asymptotic 
formula [H, 10.4.59] that 

Aiix + iy)- 1 = 0(e- cy3/2 ) (3.3) 

for some c > 0, uniformly for |x| < A for any fixed A and \y\ > 1, say 
(to avoid the zeros of Ai). By differentiation (Cauchy's estimates, see e.g. 
[la . Theorem 10.25]), it follows that the same holds for all derivatives of 
Ai(x + iy) _1 , and thus all derivatives o£g(t) decrease rapidly. In particular, 
g belongs to the Schwartz class S of rapidly decreasing functions on R; since 
this class is preserved by the Fourier transform, also g G S (see e.g. fl7l . 
Theorem 7.7]). In particular, g is integrable. 
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The characteristic function of V is the Fourier transform f(t), and thus by 
(|3.ip . (|3.2|) and standard Fourier analysis (see e.g. [13, Theorems 7.7-7.8], 
but note the different normalization chosen there), with g(t) := g(—t), 

= /(*) = g^C*) = 2 ' 2tt^*^ = 4^*^ 

i r°° ds i r io ° dz 

/ <Ls u ~ (3.4) 



27T J.oo Ai(i(t + s))Ai(is) 2iri J_ ioo Ai(z + ii)Ai(z) ' 

This is also given in Groeneboom (l3l . Lemma 2.1] (although with typos in 
the formula). The last integral is taken along the imaginary axis, but we 
can change the path of integration, as long as it passes to the right of the 
zeros a n of the Airy function (which are real and negative), for example a 
line Rez = b with b > a\ = — \a\\. 
We pause to note the following. 

Theorem 3.1. The moment generating function Ke tv is an entire function 
of t, and is given by, for any complex t, 

1 /-6+ioo j 

Ee ty = - — \ \ / \ i (3-5) 

2mJ b _ ioo Ai(z + t)Ai(z)' V ' 

for any real b with b > a\ and b + Re t > a% . 



Proof. The density function f(t) = exp(— 1 3 /3) as t — > ±oo by 11, Corollary 
3.4(iii)]. Hence, Ke tv < oo for every real t, which implies that Ke tv is an 
entire function of t. The formula (|3.5p now follows from (|3.4p by analytic 
continuation (for each fixed b). □ 

By differentiation of (|3.4p (or (|3.5p ) we obtain, for any n > 0, 



E 



_ , x „ d n . , i r 00 d n / i \ dz . , 



By integration by parts, this is generalized to: 
Theorem 3.2. For any j, k > 0, 



j+k _ r & ( 1 



2^i J_ ioo dzJ VAi(z)y dz fc VAi(z) ' 



Proof. For j = 0, this is (13.6p . 

If we denote the integral on the right-hand side of (13.7P by J(j, A;), then, 
for j, > 0, integration by parts yields A(J, k) = —A(j — 1, k + 1), and the 
result follows by induction. □ 

Since KV^ +k = EV k+ i, we see again by symmetry that KV n = when 
n is odd. For even n, it is natural to take j = k = n/2 in (|3.7p . For small 
n, this yields the following examples. First, n = j = k = yields 
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as noted by Daniels and Skyrme this is easily verified directly, since 
7rBi(z)/Ai(z) is a primitive function of 1/Ai 2 , see [2j. 
Next, for n = 2 and n = 4 we get 

2 -1 f ico / d / 1 \V -1 f 00 Ai'(z) 2 , . x 

EF 2 = / ( ~ i — ( "J — 7 — -r I I d z = / a 7 L dz (3.9 

27TiJ_ ioc \dz\Ai(z)J J 2vri 7_ ioo Ai(z)4 v ' 

and 

2vri 7_ ioo V dz 2 VAi^yy 

1 f ico f z 2Ai'(z) 2 \ 2 , . 

= Si/_(-AiW + w) ^ <3 - 10) 

Remark 3.3. Since Ai (z) = zAi(z), it follows by induction that the m:th 
derivative -^r(Ai(z) _1 ) can be expressed as a linear combination (with 
integer coefficients) of terms 

(3 ' u> 

with j, k > 0, 2j + < m and £ = k + 1. 

4. Some combinatorics of Airy integrals 

Inspired by Remark 13.31 we define in general, for any integers j,k,£ > 
with £ > k, 

1 f io ° ziA\'(z) k 

(The integrand decreases rapidly as z — > ±ioo because £ > k, by (|3.3|) and 
Ai'(^)/Ai(z) ~ —z 1 ! 2 as |z| — )• oo in any sector | argz| < tt — e with e > 
[l], 10.4.59 and 10.4.61]; thus the integral is absolutely convergent.) Then, 
recalling Ai"(z) = zAi(z), for any j, k > and £ > k, 

1 r ico d ^'Ai'(z) fe dz 

= iJ(j - 1, fc, + kl(j + l,k-l,£-l)- £I(j, k + l,£+l), (4.2) 

where we for convenience define I(j,k,£) = for j < or k < 0. Conse- 
quently, for j,k>0 and £> k, 

k + 1,1 + 1) = J jl(j - 1, fc, *) + + l,k — l,£—l), (4.3) 
or, for j > 0, and £ > k > 1, 

'^/(i- W-i), fe = i, 

- 1, A; - 1,£ - 1) + ^Efl(j + l,k — 2,£—2), k>2. 



I(j,k,£) 



(4.4) 
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By repeatedly using this relation, any k,£) may be expressed as a ratio- 
nal combination of I(p, 0, 1 — k) with < p < j + k/2. For example, 

1(0,1,3) =0; (4.5) 

7(0,2,4) = |/(1, 0,2); (4.6) 

1(1, 2, 4) = i/(0, 1, 3) + i/(2, 0, 2) = |/(2, 0, 2); (4.7) 

7(0, 4, 6) = §7(1,2,4) = ±7(2,0,2). (4.8) 

5. Back to moments 

Proof of Theorem By combining (|3.6p and Remark 13.31 we can express 
any moment E V n as a linear combination with integer coefficients of terms 
k,k + 2) with 2j + k < n. By repeated use of (14. 4p (see the comment 
after it), this can be further developed into a linear combination with rational 
coefficients of terms 0, 2) with < j < n/2. 
For example, by (j3.9|) and ()4.6p . 



E 



1 1 f iao z 

= -'(0.2.4) = --7(1,0,2) = --J_^ — iz , (5.1) 

which yields (fl3j) by (|2T2"|) . 

To continue with higher moments we have next, by (I3.10P and (|4.7p ~ (|4.8p . 

EK 4 = 1(2, 0, 2) - 4/(1, 2, 4) + 4/(0, 4, 6) 

l-| + |)/(2,0,2) = 1/(2,0,2) 

dz. (5.2) 



7 Z" 100 z 2 



30m J_ ioo Ai(z) 2 

Similarly (using Maple), 

1 r°° (26-31z 3 )/21 
27ri i_ ioo Ai(z 

In general this procedure yields, for some rational numbers 



E V 6 = / V ^. (5.3) 



Er = EM(i0.2)^/_ itj |fd, (5.4) 



where p n {z) ■= Y^j=o bnjZ^ is a polynomial of degree at most n/2. By (12. 2p . 
this is equivalent to the more general (jl.4[) . □ 

For odd n, we already know that EV n = 0, so we are mainly interested 
in p n for even n. The polynomials Po,P2,Pi,Pe are implicit in (13. 8D . (15. 1|) . 
(|5.2p and (|5.3[) , and some further cases (computed with Maple) are given in 
TableU 
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Po(z) = 1 

P2(z) = ~\z 

p 6 (z) =-§, 3 + i 
ps(z) = ^ - ^ Z 
Plo{z ) = -^z 5 + ^z 2 

„ ( v \ _ 1414477 J> _ 2419532 3 , 198947 
PY1\Z) — 1365 & 273 * " r 1365 



Table 1. The polynomials p n (z) for small even n. 



Remark 5.1. We can see from Table Q] that (for these n) p n {z) contains 
only terms z 3 where j = n/2 (mod 3). This is easily verified for all even n: 
a closer look at the induction in Remark 13.31 shows that only terms (|3,lip 
with 2j + k = m (mod 3) appears, and the reduction in (j4.4|) preserves 2j+k 
(mod 3). 



6. Problems and conjectures 

The proof above yields an algorithm for computing the polynomials p n (z), 
but no simple formula for them. We thus ask the following. 

Problem 6.1. Is there an explicit formula for the coefficients b n j, and thus 
for the polynomials p n (z)7 Perhaps a recursion formula? 

We have computed p n {z) for 1 < n < 100 by Maple, and based on the 
results (see also Table [T]), we make the following conjectures. 

Conjectures 6.2. (i) p n (z) = for every odd n. 

(ii) p n {z) has degree exactly n/2; i.e., the coefficient b n ^ n /2 of z n l 2 is 
non-zero. 

(iii) These leading coefficients have exponential generating function 

OO T 

n\ sinh x 

n=0 



Of these conjectures, (i) is natural, since we know that KV n = for odd 
n, and (ii) is not surprising. The precise conjecture (|6.ip is perhaps more 
surprising. We have verified that the coefficients up to a; 100 agree, but we 
have no general proof. 

The simple form of (|6.ip suggests also the following open problem. 

Problem 6.3. Is there an explicit formula for the generating function 

En=0Pn(z)x n ? 



8 



SVANTE JANSON 



References 

[1] M. Abramowitz & I. A. Stegun, eds., Handbook of Mathematical Func- 
tions. Dover, New York, 1972. 

[2] J. R. Albright & E. P. Gavathas, Integrals involving Airy functions. 
Comment on: "Evaluation of an integral involving Airy functions" by 
L. T. Wille and J. Vennik. J. Phys. A 19 (1986), no. 13, 2663-2665. 

[3] D. Anevski and P. Soulier, Monotone spectral density estimation. Ann. 
Statist. 39 (2011), no. 1, 418-438. 

[4] A. D. Barbour, A note on the maximum size of a closed epidemic. J. 
Roy. Statist. Soc. Ser. B 37 (1975), no. 3, 459-460. 

[5] A. D. Barbour, Brownian motion and a sharply curved boundary. Adv. 
Appl. Probab. 13 (1981), no. 4, 736-750. 

[6] H. Chernoff, Estimation of the mode, Ann. Inst. Statist. Math. 16 
(1964), 31-41. 

[7] H. E. Daniels, The maximum size of a closed epidemic. Adv. Appl. 

Probab. 6 (1974), 607-621. 
[8] H. E. Daniels, The maximum of a Gaussian process whose mean path 

has a maximum, with an application to the strength of bundles of fibres. 

Adv. Appl. Probab. 21 (1989), no. 2, 315-333. 
[9] H. E. Daniels & T. H. R. Skyrme, The maximum of a random walk 

whose mean path has a maximum. Adv. Appl. Probab. 17 (1985), no. 1, 

85-99. 

[10] P. Groeneboom, Estimating a monotone density. Proceedings of the 
Berkeley conference in honor of Jerzy Neyman and Jack Kiefer, Vol. II 
(Berkeley, Calif, 1983), Wadsworth, Belmont, CA, 1985, pp. 539-555. 

[11] P. Groeneboom, Brownian motion with a parabolic drift and Airy func- 
tions. Probab. Theory Related Fields 81 (1989), no. 1, 79-109. 

[12] P. Groeneboom, The maximum of Brownian motion minus a parabola. 
Electronic J. Probab. 15 (2010), paper 62, pp. 1930-1937. 

[13] P. Groeneboom, Vertices of the least concave majorant of Brownian 
motion with parabolic drift. Electronic J. Probab. 16 (2011), paper 84, 
2234-2258. 

[14] P. Groeneboom & J. A. Wellner, Computing Chernoff's distribution. J. 

Comput. Graph. Statist. 10 (2001), no. 2, 388-400. 
[15] S. Janson, G. Louchard Sz A. Martin-L6f, The maximum of Brownian 

motion with parabolic drift. Electronic J. Probab. 15 (2010), paper 61, 

pp. 1893-1929. 

[16] W. Rudin, Real and Complex Analysis. McGraw-Hill, London, 1970 
[17] W. Rudin, Functional Analysis. 2nd ed., McGraw-Hill, New York, 1991. 

Department of Mathematics, Uppsala University, PO Box 480, SE-751 06 
Uppsala, Sweden 

E-mail address: svante . janson@math.uu. se 
URL: http : //www . math . mi . se/ ~svante/ 



